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1. Introduction 

How to prepare a quantum state? It is a nontrivial and important problem to be 
tackled. In fact, various interesting and peculiar phenomena are predicted on the basis 
of highly nonclassical states, and entanglement plays key roles in quantum information 
protocols [T] . They all rely on the generation of nontrivial states and are not realized 
without establishing the strategies for the preparation of such quantum states. 

Generally speaking, we try to drive a quantum system to a specific state by a series 
of operations, e.g., applications of external fields to transform its state, measurements 
to project it onto a particular configuration, and so on. A generic mechanism was 
found to extract a pure quantum state from a given arbitrary (mixed, in general) 
state, by simply repeating the same measurement on an ancilla system interacting 
with the target system [2j[3]. Such a mechanism is interesting in itself and is even 
indispensable when direct operations on the target quantum systems are not allowed 
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or unavailable. The repeated measurements on the ancilla can be regarded as an 
indirect (POVM) measurement on the target system which, under proper conditions, 
allows us to drive the latter toward the desired pure state. This idea was applied 
to the initialization of qubits [3] , extraction of entanglement [3113] and a nonclassical 
state [5], and to establish entanglement between separated qubits [MS]- 

In those schemes, a pure quantum state is obtained from an arbitrary initial 
configuration, only when the ancilla system is repeatedly confirmed to be in a specific 
state by all the measurements performed during the protocol. That is, they are 
probabilistic schemes characterized by a success probability strictly less than 1. The 
primary motivation of the present work is to persue a scheme which would allow one 
to reach the target state with probability 1, or at least, with probability arbitrarily 
close to f, independently of the initial conditions. To achieve such a goal, we take 
inspiration from an approach recently introduced in Ref. |10j , in which a target system 
is indirectly controlled by making it interact with a sequence of (properly prepared) 
ancillas, which are then discarded. As in the case of Ref. [TU], our finding relies on 
a useful property of quantum channels. Namely, we make use of the fact that under 
proper conditions (see, for instance, Refs. [HJ[12] and references therein) repetitive 
applications of the same map drive the system toward a fixed point, independently of 
its initial configuration (mixing property). 

As a nontrivial example of such scheme, we concentrate on a prototypical setup 
which has been extensively investigated in the literature recently [7H91 IT3HT5] . Here, 
two non-interacting target qubits A and B sit at fixed distance from each other along 
a ID channel as sketched in Fig. [TJ The goal of the scheme is to drive A and B into an 
entangled state with the help of a (flying) qubit, which is sent through the ID channel 
and is detected after it has been scattered by the targets. In the simplest configuration 
considered so far, the latters are supposed to be initially in a (known) separable state, 
while the ancilla qubit is prepared in an appropriate spin state before being injected 
into the setup. Under these assumptions, it has been shown that entanglement between 
A and B can be generated in a probabilistic fashion by a simple post-selection of the 
spin state of the scattered ancilla . Schemes which do not require preparation 

of the initial state of the target qubits to a specific state were also proposed, in 
which entanglement is extracted after repetition of scattering + post-selection [THS] . 
Such improved protocols however are still probabilistic, as they produce the desired 
entanglement only with a certain success probability. In contrast, the approach we 
present here generates a maximally entangled state between A and B from an arbitrary 
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Figure 1. A schematic drawing of the setup. An ancilla qubit X is sent to two 
fixed qubits A and B, with a certain wave vector k. After being scattered by the 
delta-shaped potentials produced by A and B, we check whether it is reflected 
or transmitted. Neither the preparation of the spin of the incident X nor the 
spin-resolved detection of the scattered X is required. 
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initial state with probability 1. Furthermore, it requires neither the preparation nor 
the post-selection of the spin state of the ancilla qubits. 

The paper is organized as follows. In Sec. [3J we introduce the setup and give 
some preliminary definitions. Section [3] introduces the protocol and states the main 
result of our work. The proof of the latter is then provided in Sec. [4] Section [5] is 
devoted to analyze the robustness of the scheme, while conclusions and remarks are 
given in Sec. |6] 



2. Setup 

Our setup is sketched in Fig. [1] Two qubits A and B are fixed at x = —d/2 and d/2, 
respectively, along a ID channel. They do not interact directly with each other, while 
we wish to establish entanglement between them. To do so, we send flying ancilla 
qubits X as "mediators" and let them scatter with A and B. As in Rcfs. [7H MT3"I[T4] . 
we assume the system to be described by the following Hamiltonian: 

+ g{v™ . -d/2), (2.1) 

where x and p are the position and the momentum of X in ID, the operators cr^ 
(J = X, A, B) represent the Pauli operators of the spins, and the potentials produced 
by A and B are represented by the delta-shaped potentials. According to Eq. (12. ip . 
the spin of X interacts separately with A and B through the Heisenberg-type coupling 
during the scattering. This Hamiltonian has been proposed to effectively model the 
coupling between electrons occupying the lowest sub-band and magnetic impurities 
placed along a quasi one-dimensional (ID) wire, such as a semiconductor quantum 
wire |16j or a single- wall carbon nanotube |17j . where electrons flow. 

The particle X is sent from the left with a fixed incident wave vector k > and 
scattered by A and B. The matrix elements of the scattering operator S are given 
by (HUH] 

(k'C\S\K) = e- ikd [S(k' - k)(('\T k \C) + 6(k> + k){C\R k \0], (2.2) 

where \k) is the eigenstate of the momentum operator p of X belonging to its eigenvalue 
hk, and |0 represents a spin state of XAB. The operators Tk and Rk describe the 
changes provoked in the spin state of XAB when X is transmitted to the right and 
reflected to the left, respectively. They satisfy the unitarity condition 

TlT k + RlR k = l XAB , (2.3) 

and are given by 

T k = e M [a fc (l - 4in k )P_ + (a k Qi + (3 k Qs)P + 

- a fc fi|(l - e 2lM )(P_ - '3QiP+ -K+ + A'_)] , (2.4a) 

R k = T k e~ ikd - 1 - ifi fc (l - e 2M ) S^a k n 2 k (l - e 2lkd )P_ + (2a k Qi - /3 k Qs_)P + 

+ ia fc [l + 30^(1 - e 24M ) - 4iQ k P + ](K + + A„)|, (2.46) 
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(1 - 4ifi fc ) + 2fi£(l - 6iCt k j{l - e Zlka ) + 9fi£(l - e 2 



ak n an s . nn2/, c „r> \fi „2ikd\ , m 4h „2ifcd\2' (2.5a) 



^ ~ (1 + 2zr! fe ) - n|(l - e 2ifcd ) ' ° fc ~ t£' (2 ' 5&) 
In the above expressions, 

arc the projection operators on the singlet and triplet sectors of A and B, respectively, 
while 

Ql = l P+ + \*W.{ fr W+oM), (2.7a) 

2 6 

Q,=P- + ip + -^«> (A) +<^) (2.76) 

are those on the spin-| and spin-i sectors of XAB, respectively [TS] (note that they 
are all commuting with each other and satisfy QaP_ = P_Qa = 0)0 The other 
operators 

K^^I.Ef, (2.8) 

defined with 

S CA*> = 1 [(<T (A) _ ^(B)) + l{(T {A) x fj.{B) 

= P + (<t^ - <r( B >) = P+i(o-< A > x o-W) 

= (<T< A ) - CT ( B ))P_ = i(<7^ X CT ( B ))P_ 

= E^ S)t , (2.9) 

are responsible for the transitions between the singlet and triplet sectors of A and B, 
with the only nonzero elements P±Z [ ± B) P T ^ 0. 

It is pointed out in Ref. [14] that this system exhibits interesting resonant 
transmissions controlled by the entanglement in A and B: for instance, when A and 
B are in the singlet state \^~)ab = (ITDab — |4-t) ab) I 'V%, the potentials produced 
by A and B look "transparent" for X sent with a resonant wave vector k n satisfying 
the resonance condition k n d = rnr (n = 1,2,.. .). Sec Fig. [2] More explicitly, the 
transmission and reflection operators T k and Rk at the resonance points are given by 



Tk n = (-1)" 



Qi Qs. 

p I | 2 | 2 1 p 

' - 4ifl k 1 + 2iQ k 



(2.10a) 



Rk " = (i^ * - TT^t g i) P +' (2 ' 10&) 

showing that X is perfectly transmitted without spin flip when A and B are in the 
singlet state \^/~)ab- 

J Throughout this paper, the unit operators are often omitted as lx CE5-P± — > P± > lx ® f^ 1 *' ® 1b — > 
<t< a ), 31 ab -> 3, etc. 
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Figure 2. Transmission probability of the ancilla qubit X when it is injected 
from the left with its spin prepared in |t)x- The plot shows the probability of 
detecting the ancilla on the right detector after the scattering by A and B. Here A 
and B arc assumed to be initially in the singlet state \^~)aBi so that the quantity 
plotted is nothing but Tr{T fc (|t)x(tl ® \^~)ab{^~ \) T t}- Its dependence on the 
incident wave vector k of X and the coupling constant g is shown. It exhibits 
resonances at k = nn/d (n = 1, 2, . . .) for any g ^ 0. 



3. Protocol 

To construct our scheme, we make use of the resonance condition detailed in the 
previous section. Similar approaches have been explored in Rcfs. [7 9 . Specifically, 
in Rcf. [7] , Ciccarello et at exploited repetition of the scattering of X with a resonant 
momentum followed by an appropriate post-selection on X, to extract the singlet state 
I^^ab from AB with a probability which depends upon the initial state of the system. 
In contrast, in Ref. [HUH], a scheme is proposed to extract the singlet state from an 
arbitrarily given initial state of AB, in which neither the preparation of the spin state 
of X nor its post-selection is required. However, it is still a probabilistic scheme, since 
the singlet state is extracted by the successive post-selections of transmitted events. 
The protocol presented here solves this problem, allowing one to produce a maximally 
entangled state of A and B from their arbitrary initial state with probability 1. The 
scheme remains free from the preparation and the post-selection of the spin state of 
X. 

Here comes our protocol (see Fig. 

0. The initial state of A and B is arbitrary and is in general a mixed state p. 

1. We send X with its spin arbitrary from the left to A and B with a resonant 
wave vector k n = mr/d (n — 1, 2, . . .), and see if it is transmitted to the right or 
reflected to the left, irrespectively of its spin state. 

2. If X is detected on the right (transmitted), we proceed to the next round (to step 
1). 

3. Otherwise (reflected), we send another X from the left with its spin randomly 
chosen with an off-resonant wave vector q (perfectly polarized incident spin of X 
is not recommended). We do not check anything after this scattering; we just 
proceed to the next round (to step 1). 
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Figure 3. Flowchart of the protocol. 



4. We repeat this routine (steps 1 to 3) many times and end up with the singlet 
state \^~)ab in A and B with probability 1. 

This is a feedback approach, where we take different actions depending on the 
outcome of the resonant scattering at step 1. Notice however that, for implementing 
the feedback, no additional element is required: we have only to set the incident wave 
vector of the ancilla off-resonant. 

It is also worth stressing that, at step 1, the choice of the incident spin of X 
is irrelevant and we can just choose it arbitrarily. At step 3, on the other hand, we 
can choose it randomly, but if the incident spin of X is perfectly polarized always in 
the same direction at every cycle, the scheme does not work. In order to make the 
following analysis simpler and transparent, however, we take unpolarized spin for both 
steps, represented by the completely mixed state lx/2. 

A single cycle (steps 1 to 3) changes the spin state p of A and B in the following 

way: 

p^Mp = (T kn +S q 1l kn )p, (3.1) 

where 

%p = Tr x {T k (t x /2 ® p)T*}, (3.2a) 
K k p = Tr x {R k (l x /2®p)Rl}, (3.26) 

and 

S k =T k +K k . (3.3) 

On the basis of the unitarity (|2.3p . the map M. is trace preserving, 

Tr AB {Mp} = 1, (3.4) 

meaning that all possible outcomes are properly kept at each step and there is no 
selection of specific detection events. 

We will see that the repetition of the above cycle leads A and B into the singlet 
state, 

M N P ^ (JV ->•«>), (3-5) 

irrespectively of their initial state p, where N is the number of cycles of the protocol. 
In other words, we are going to prove that M. is "mixing" [TT] with its fixed point 
given by |*-)ab(*~|- 
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4. Proof 

In order to verify the claim (|3.5j) . we invoke an important result on mixing channels, 
which states that a CPT (completely positive and trace-preserving) map M. is mixing 
if it has a unique fixed point that is a pure state, e.g., see Ref. [TT]. We remind that 
the fixed points of M. are defined as those input states p* which are left invariant by 
the action of the channel, i.e., 

Mp*=p*. (4.1) 

In our case, it is immediate to check that \^~)ab(^~ | is a fixed point of Ai: indeed, 
as is clear from Eqs. (|2.10a[) and (|2.106[) . \^~)ab(^~ | is preserved by 7fe„ while 7Zk n 
yields nothing (reflection does not occur with \^^)ab(^^\)- Thus, since \ 1 $>~)ab( 1 &~ \ 
is a pure state, it follows that the only thing we need to verify in order to prove the 
mixing (|3.5[) is that there exists no other fixed point of the channel M. . 

Assume then that there exists another fixed point p*, that is different from 
\^~)ab(^~ I- By definition, it must satisfy the identity 

M P * = {%„ + S q K kn )p. = p* . (4.2) 

This expression can be simplified by noticing that the following identity holds at 
resonances: 

P-T kn =V-, (4.3) 

where V± are the superoperators associated with the projections onto the triplet and 
singlet subspaces, respectively, 

V±p = P± P P ± . (4.4) 

Indeed, taking Eq. (|4.3[) into account, a necessary condition for Eq. (|4.2[) is given by 

V-S q TZ kn p, = 0. (4.5) 

Look at Eq. (|2.106[) again: it shows that TZk n cuts out the singlet sector and acts 
on the triplet components. Since we are assuming that p* is different from the 
singlet state \^>~)ab(^~\, we surely find some component lZk n p* ^ in the triplet 
sector. Therefore, if S q is such a map that certainly couples any triplet components 
to the singlet sector, the condition (|4.5[) is never satisfied except for the singlet state 
I^^ab^ - !- Consequently, by contradiction, the singlet state is proved to be the 
unique fixed point of the map M.. 

The condition for such a map S q is expressed as 

TTA B {V-S q V+p} > 0, (4.6) 

for any state p ^ \'$>~)ab( 1> &~ |- By inserting the explicit expressions of T q and R q 
given in Eqs. (|2.4al) and (|2.4fe|) . it reads 

Tv AB {VS q V + p} = W q Tr AB {P+p}, (4.7) 

with a coefficient 

4n 2 q + |1 - 2iQ q + 30^(1 - e 2?;<jd )| 2 j , (4.8) 



W q = \a q \ 2 n 2 q \l ~ e 2i « d \ 2 



which is nonvanishing for any q =/= rm/d (n = 0,1,2,...) and for any state p =/= 
\^~)ab(^~\- That is, S q with an off-resonant wave vector q is a map that surely 
couples any triplet component of A and B to the singlet sector and excludes the 
existence of p* satisfying Eq. (|4.5[) and different from the singlet state \^^) ab(^~\- 
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See Fig. |4j where the coefficient W q is plotted as a function of q and g. They are 
actually nonzero except at the resonant wave vectors q — nir/d (n= 1,2,.. .). 

The reason why A and B are attracted into the singlet state by the repeated 
applications of M. is the following. The singlet state \^>~)ab(^~ | is the fixed point 
of the map Ai, and the singlet component in the state of A and B remains there. As 
for the triplet components, they provoke the reflection of X with a certain probability 
at step 1. Once X is found to be reflected, the qubits A and B are "shuffled" by the 
subsequent off-resonant scattering S q , which creates a singlet component. In total, 
the probability to find A and B in the singlet state is increased by the single cycle. 
Outflow of the probability from the singlet sector is absent, while inflow is present. 
This feature leads the system into the singlet state \^~)ab(^~ |H 

The scheme surely works. As a direct check, the average fidelity F(N) of the 
protocol is shown in Figs. [5] and [B] It is obtained by computing the fidelity between 
the generated state Ai N p and the target \^~)ab('^~\, which is averaged over all 
possible choices of the input state p, i.e., 

F(N) = Ty ab {V-M n P } = Ti ab {V-M n (1ab/4)}, (4.9) 

where • • • stands for the average over the input state p and 1,4s is the identity 
operator on AB. From these plots, it is evident that as the number N of protocol 
cycles increases, the average fidelity F(N) approaches asymptotically to 1, as long as 
the wave vector q is properly set off- resonant. This implies that for an average choice 
of the input state p, the state M N p approaches \^~)ab{^~ |- By linearity, this also 
implies that the same result should hold for each input 

p(U Therefore, Figs. and |6] 

§ One may wonder that a reflection event at step 1 suddenly projects out the singlet component and 
keeps A and B from approaching the singlet state. As the cycle is repeated, however, the average 
probability to find A and B in the singlet state increases (as explained above), and the average 
probability for such a reflection to occur is accordingly reduced: the chance of the loss of the singlet 
component becomes less and less likely, and the singlet component keeps on growing on average. 
A simple statistical argument can then be used to claim that, since on average these events are 
suppressed, the probability of generating a specific trajectory which does not have this property is 
also suppressed. 

|| As a matter of fact, one can easily verify that in case of pure fixed point, asymptotically optimal 
average fidelity is a sufficient condition for the mixing property of a CPT map. 
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Figure 5. Average fidelity F(N) defined in Eq. | |4,9| I as a function of the number 
N of protocol cycles. The five curves refer to different choices of the resonant wave 
vectors, namely k„d/ir = f , 2, 3, 4, 5 from top to bottom. The other parameters 
are qd/n = 2.5 and g = h 2 Tr/md. 




Figure 6. Average fidelity F(N) defined in Eq. I I4.9I I. (a) Dependence on q with 
k n = Tt/d and g = % 2 ir/md. (b) Dependence on g with k n = n/d and q = 2.5n/d. 



provide an alternative proof of the mixing property of A4. 

The speed of the convergence is controlled by the following two factors: (i) the 
reflection of X sent with a resonant wave vector k n = rnr/d (n = 1, 2, . . .), at step 1, 
due to the presence of triplet components in AB, and (ii) the transition of AB from the 
triplet sector to the singlet by the off-resonant scattering with q ^ mr/d (n = 1, 2, . . .), 
at step 3. The resonant scattering does not bring any triplet components of AB to 
the singlet sector at step 1. However, the reflection of the resonant X triggers us to 
proceed to the off-resonant scattering at step 3, which brings some triplet components 
of AB to the singlet sector. Therefore, the larger is the reflection probability of the 
resonant X and the stronger is the transition from the triplet sector to the singlet 
by the off-resonant scattering, the faster is the convergence to the fixed point. The 
former is controlled by Rk n in Eq. (|2.106[) and the latter by W q in Eq. (|4.8|) . In 
particular, the reflection probability is smaller for a higher incident wave vector k n . 
See Slk n in the numerator of Rk n in Eq. (| 2 . 1 &[) . W q in Eq. (14. 8} is also proportional 
to Vl 2 q and is a decreasing function of q (for not too small q) (apart from the oscillation 
due to resonance). The convergence is thus slower with higher incident momenta, as 
demonstrated in Fig. [5] 
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5. Robustness 

In this section, we analyze the robustness of the proposed scheme. We start by 
considering errors in the preparation of the ancilla qubits. Later, we analyze the 
effect of inefficient detectors. 

5.1. Errors in the incident momentum 

A key ingredient of our entanglement protocol is the resonant tunneling condition of 
the flying ancillas we impose at step 1. To see what happens if one fails to enforce such 
constraint, we consider two alternative scenarios. First, we analyze the case in which 
the source producing the ancilla qubits X is affected by a systematic error that forces 
it to produce a monochromatic sequence of particles which enters the ID channel with 
a constant wave vector k, which is not resonant (error by deviation from resonance). 
Then, we consider the situation in which the same source is affected by fluctuations 
which prevent it from producing monochromatic signals (error by fluctuation of the 
incident momenta). For both scenarios, we compute the overlap between the final state 
of AB after N protocol cycles, and the target singlet state. As one might expect, the 
case of the systematic error is much more detrimental for the performances of the 
protocol, with average fidelities which drops below 50% already for small deviations of 
the impinging momenta. On the other hand, the scheme appears to be more resilient 
to fluctuation errors. 

5.1.1. Deviation from a resonance point Assume that, at step 1 of the protocol, the 
ancilla qubits X enter the ID channel with fixed wave vector k, which is not necessarily 
at resonance. Following the derivation in Sec.[3j one can easily verify that, after each 
protocol cycle, the transformation of the state of AB can still be described as in Eq. 
(|3.ip . but with the superopcrator Ai replaced by the CPT map 

M k , q = T k +S q TZ k . (5.1) 

Ideally, k should be set at a resonance k n = nir/d (n = 1,2,...), while q should 
be off-resonant. As shown in the previous section, these assumptions are sufficient 
to guarantee that the map M = Mk n .q is mixing with the singlet state as its fixed 
point. For k ^ k n , however, this is not necessarily true, posing the problem on how 
to compute the state of AB in the asymptotic limit of large N (if M.k, q is not mixing, 
limjv->oo might not be well defined with the system continuously oscillating 

between different configurations |llj). For the sake of simplicity, however, we will 
neglect this issue in the following, assuming the mixing property to hold anyway. Even 
though we do not have a formal proof of this property, such assumption is strongly 
supported by a series of numerical and theoretical evidences. We remind in fact 
that the mixing property of a CPT map is ultimately related with its spectrum |llj . 
Specifically, a necessary and sufficient condition for mixing is the existence of a finite 
gap between the largest and second largest eigenvalues of the map. As shown in Fig. 
this seems to be the case for M.k, q - Furthermore, since Mk, q depends continuously 
on k, its spectrum is a continuous function of fc, and hence, so is its mixing property. 
That is, there is at least a neighborhood of k n = nir/d (n = 1,2,.. .) such that all 
k G (fc„ — e, k n + e) give rise to mixing maps. Finally, it is known that the set of 
non-mixing channels form a subset of zero-measure in the set of CPT maps, so it is 
highly unlikely to have M.k, q non-mixing [TP] , 



Entanglement by repeated on- and off-resonant scattering of ancilla qubits 



11 



1.0 



0.8 
— 0.6 





^0 t\ 
n i i 
J ii ' i 












II ' 

f 








= 2.57T 
g = Ti'ix jmd 



1.0 — ' — ' — ' — ' — ' — i — i — i — ■ — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i 
1 2 3 4 5 

kd/lT 



Figure 7. The magnitudes of the largest and the second largest eigenvalues, Ao 
and Ai, of the map M k ,q as a function of k with q = 2.5n/d and g = h 2 ir/md. 
The presence of a gap between them is a necessary and sufficient condition for 
mixing |11| . 



With the above considerations in mind, we identify the state of AB after TV 3> 1 
protocol cycles for generic k with the fixed point p*(k) of the map M k ,q, i-e., 

M k , qP *{k) = p*{k). (5.2) 

Interestingly enough, even without solving Eq. (|5.2[) explicitly, it is possible to derive 
a concise formula for its fidelity with respect to the singlet state, 

F t (k)= Tr AB {V- p^k)}. (5.3) 

To see this, we first notice that the transitions between the singlet and triplet sectors 
of A and B induced by a single scattering event are described by the following 2x2 
matrices: when X is transmitted, 



with 



Tr AB {V-%p} \ = / T k — T h 7+ \ ( Tr AB {V-p} 
Tv AB {V+%p} )~{ T+- T++ ) { Tr AB {V +P } 



and similarly, when X is reflected, with 
'ft 



K~ + = = \a k \ 2 n 2 k \l ~ e 2tkd \ 2 \l - 2iVL k + 3Q|(1 - e 2ikd )\ 2 , (5.6b) 



(5.4) 



r k 7 + = -T+- = 4|a*| 2 fi£|l - e 2ikd \ 2 , (5.56) 



T k — = \a k \ 2 \l - 4itl k - - e'T, (5.5a) 
1, 
3 

T++ = \\{a k + 2 fa) + 3a k Q 2 (l - e 2ikd )\ 2 

+ \\{a k - fa) + 3a k Sl 2 (l - e 2ikd )\ 2 , (5.5c) 



K-- = \l- a k (l - 4iO fc ) + a k n 2 (l - e 2M ) + 6ia k n 3 k (l - e 2ikd ) 2 \ 2 , (5.6a) 



K++ = -|3 - K + 2/3 fe ) + ifi fc [2(a fc - &) + 3ia k n k ](l - e 2ikd )\ 2 

+ \\{oc k ~ fa) - in k [(2a k + fa) + 3ta k n k ](l - e 2lkd )\ 2 . (5.6c) 
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In this notation, W q denned in Eq. (|4.8j) is expressed as W q = T q + +lZ q + , and the 
2x2 matrix for S q reads 

/ 1 - 3W q W q \ 
S q =( . (5.7) 

V ZW q 1-W g J 

[By abuse of notation, we use the same symbols for the corresponding 2x2 matrices, 
e.g., S q for the 2x2 matrix in Eq. (|5 . T[) .] Combining these expressions, the matrix for 
A4k,q is also constructed according to Eq. (|5.1|) . In particular, at resonances k n = nn/d 
(n = 1,2,...), we have 




with 

[1 + 160|J(1 + 4Q^ 

and for the ideal map, 



/ 1 W q V n \ 
M=M kn . q =\ . 5.10 

V o i - vK 9 y„ y 

Now, by noting Eq. (|5.2[) . "P+ = 1 — "P_, and Tr^B /)»(^) = 1; the definition of the 
fidelity ()5.3|) is arranged as 

F„(fc) =Tr AB {P_^(fc)} 

= Tr AB {7>-.M M p*(fc)} 

= -Mfc,, Tr AB {P-p*(k)} + Mfc + Tr AB {P + ^(fc)} 

= + (Mfc " - Mi" +) TWP_p*(A0} 

= Mfc+ + (Mfc--^+)F,(fc). (5.11) 

Therefore, as long as 1 — + > (which is assured by > or 

M~^~ q < 1), one gets a concise formula for the fidelity, 

*<*>-, -xgLa- ( " 2) 

In Fig. [51 we report its plot as a function of the incident wave vector k for a fixed g: as 
anticipated, the fidelity F*(k) drops below 50% as fc deviates from a resonance point 
k n = nir/d (n= 1,2,...) by less than a percent. 



5.1.2. Wave packet or fluctuation of the incident momenta Consider now the case in 
which the source emitting X injects a stream of non-monochromatic ancillas into the 
ID channel. Specifically, let tp(k) and <fi(q) represent the wave packets in momentum 
space of the particles produced by the source at steps 1 and 3 of theprotocol, 
respectively. Then, the map Mk, q of the previous section is substituted bjol 

/■OO pOO 

M= dk dq\*jj(k)\ 2 \0(q)\ 2 M k , q . (5.13) 
Jo Jo 

If Wc assume that the wave packets are composed only of positive momenta. It is possible to show 
however that the presence of the negative momenta does not spoil the mixing property to be argued 
in this section. In any case, it should be reasonable to assume that such components are negligibly 
small. 
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Note that the trace over the momentum degrees of freedom is taken since the detectors 
do not resolve the momenta, and as a result, only the diagonal components with 
respect to the momenta contribute to the formula. This implies that a different type 
of fluctuation in momentum, incoherent fluctuation, is described by formally the same 
formula as Eq. (|5.13l) . Indeed suppose that the incident wave vectors k and q differ 
from run to run. Then, the state generated after N cycles of the protocol reads 

P(k uqi ),...,(k N , qN ) = M kN ,q N ■ ■ ■ M kuqi p. (5.14) 

If the fluctuations of the mementa (hi, qi) at each cycle i = 1, . . . , N are characterized 
by a common probability distribution function f(ki,qi), the state (|5 . 14[) averaged over 
the probability distribution reads 

K N ) = J (fl d h d <li f(h,qi)j P(fc ligi ),. = M N p, (5.15) 
where in this case 

/>oo />oo 

M= dk dqf{k,q)M k , q , (5.16) 



which coincides with the expression of Eq. (|5.13[) by identifying f(k, q) with 

\m\ 2 \m\ 2 - 

It is worth stressing that, differently from the case treated in the previous 
subsection, one can show that the average map A4 is mixing, provided that the 
distribution /(fc, q) overlaps with a resonant wave vector in k and has only measure 
zero in q at resonances. Indeed, split Ai into two parts as 

pk n -\-e />oo 

M= dk dqf(k,q)M k , q + M'. (5.17) 

Jk n -e JO 

Recall then that any nontrivial convex sum of a mixing map £ and something else £' 
(not necessarily mixing), 

£ = A£ + (1-A)£' (0<A<1), (5.18) 

is also a mixing map |10j . Since the first part of Eq. (|5.17[) is mixing, and has 
measure nonzero, this theorem ensures that M. is also mixing, implying that the 
mixing property of A4 is robust against the momentum fluctuation. As a consequence, 




Figure 8. Fidelity F*(fc) of the fixed point of the map A4t,q as a function of k 
for q = 2.5n/d and g = h 2 n/md. 
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Figure 9. Average fidelity F* given in Eq. I|5.19|l of the fixed point of the map 
JV[ for a Gaussian distribution of fc in Eq. I|5,20|l , (a) as a function of Afc for 
different resonant wave vectors k n d/ir = 1 (solid), 2 (dashed), 3 (dotted), with 
qd/n = 2.5 and g = h 2 ir/md, and (b) as a function of Afc and g with k„d/ir = 1 
and qd/n = 2.5. 



in the limit of infinitely many protocol cycles, the system AB is driven to the fixed 
point of the channel M. in Eq. (|5.16j) . Its fidelity with respect to the singlet state 
can now be computed similarly to Eq. (|5.12j) . yielding 

F, = Tr AB {V-M = - -~ T —, (5.19) 

1 - M — + M~+ 

where the matrix elements jCi"^ are obtained by averaging over the distribution 

f(k,q). 

In Fig. |9j the fidelity F* is plotted for a Gaussian distribution of k, centered at a 
resonance k n with width Afc, for a fixed q, 

f(k, q) (X 0(k)e- {k - k " )2/2iAk)2 S(q - q), (5.20) 

which is normalized as J °° dk dq f(k,q) = 1. As expected, the fidelity F* decreases 
from unity, as the width of the distribution Afc grows. Notably, in this case, the scheme 
seems to be quite resilient to the noise: the fidelity drops below 50% only for Afc/fc„ of 
the order of 5%. Furthermore, the decrease is slower than linear for sufficiently small 
Afc, and it is less pronounced for larger central resonant wave vector fc n and a smaller 
(but not too small) coupling constant g (although the approach to the final state is 
slower with a larger k„ , as demonstrated in Fig. ITU]) . 



5.2. Detector efficiency 

In this section, we analyze how defective detections of the scattered ancillas deteriorate 
the performances of the protocol. In particular, we consider the case in which the 
detectors in Fig. [TJ are characterized by an efficiency n < 1, i.e., they fail to report the 
arrival of a particle X with probability 1 — 77 (for the sake of simplicity, we assume 
that the two detectors have the same efficiency rj). 

To account for such events, we need to specify the action one has to take when 
the detectors fail to report the arrival of a scattered particle incident on resonance 
at step 1 of the protocol (since we do not check anything after the scattering of an 
off- resonant particle at step 3, the efficiency of the detector does not matter for this 
step). Specifically, we analyze two alternative solutions. 
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Figure 10. Average fidelity F(N) = Tr AB {V-M N (t A B/4)} as a function of 
the number N of protocol cycles, for a Gaussian distribution of k in Eq. J5 . 20t . 
The three curves refer to different choices of the resonant wave vectors, namely 
k n d/TT = 1, 2, 3, with the other parameters fixed at (Ak)d/n = 0.05, qd/n = 2.5, 
and g = H 2 n/md. 



Case I: We may simply proceed to the next round (step 1) to send the next 
particle on resonance [Fig. UlT a)]. In such a case, the map M. is modified to 

Mfi =r)M + (l-r))S kn , (5.21) 

which is still mixing by the same argument for Eq. (|5.17[) . By noting the expressions 
in Eqs. (|5.T[) and (|5.10[) . its 2 x 2 matrix which describes the transitions between the 
singlet and triplet sectors reads 

' 1 T)W q V n \ 



(5-22) 

v 1 - r)WgV n J 

and by applying the formula (|5.12[) . the fidelity of its fixed point to the target singlet 
state is shown to remain 

Fi 1) (? 7 ) = l. (5.23) 

Therefore, the protocol is still able to extract the singlet state from AB. Notice 
however that the speed of the convergence of the scheme is affected by rj < 1, as 



(a) 



(^^arbitrary state p ^ 



(b) 











resonant scattering 


No ~ 




1 -<T" rpflppt.pH? s> 

click — 
Yes 1 


off-resonant scattering | 







arbitrary state p ^ 




Figure 11. Flowcharts for Cases I (a) and II (b). 
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Figure 12. Plot of the average fidelity FW(N) = 1r AB {V-M^ )N (1as/4)} 
for inefficient detectors, with the strategy described in Case I of Sec. 15.21 N is the 
number of protocol cycles, while r) is the detector efficiency. The other parameters 
are k n = n/d, q = 2.5ir/d, and g = h 2 n/md. 



shown in Fig. [T3J This is a consequence of the fact that, in this case, the element of 
the scattering matrix (|5.22[) associated with the transition from the triplet sector to 
the singlet gets degraded: the smaller is r], the slower is the speed of convergence. 

Case II: We can proceed to step 3 (as if the ancilla X on resonance is reflected) 
to send X with the off-resonant wave vector q [Fig. [TTT b)]. In this case, the map M. 
is changed to 

=r)M + (l-r))S q S kn , (5.24) 
which is also mixing. It yields 

/i- 3( i-„h, <i-„ + ,vyw f \ 

V 3(1 - T])W q 1 - (1 - ?] + vV n )W q J 
and the fidelity of its fixed point is given by 

F * 2) W = A«~ V \V V V > ^ 
4(1 - 77) + rjV n 

which ranges between 0.25 and 1. This strategy is thus not as effective as the previous 
one in terms of fidelity. This is because the off-resonant scattering in the absence 
of guarantee that there is no singlet component in AB, which should be assured by 
the detection of a reflected particle on resonance, provokes the undesired transition 
from the singlet sector to the triplet. The presence of this leakage channel hinders the 
convergence to the singlet state. 

6. Conclusions and remarks 

We have proposed and studied a scheme for preparing a maximally entangled state 
in two non-interacting qubits, initially given in an arbitrary state. By repetition of 
resonant scattering of ancilla qubit, followed by off-resonant scattering if necessary, 
qubits A and B are driven from any initial state into the singlet state with probability 1 
(perfect efficiency). Neither the preparation nor the post-selection of the ancilla spin 
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state is required. By introducing an appropriate feedback strategy, the previously 
proposed probabilistic scheme has been turned into a reiterative scheme which leads 
the target qubits to the singlet state with probability that converges asymptotically to 
1 in the number of iterations. It is remarkable that no additional clement or technology 
is required for the feedback: we have only to set the incident momentum of the ancilla 
off a resonance point. The convergence to the unique fixed point (mixing property) 
is rigorously proved, and is shown to be robust again various types of imperfections 
in the scheme. In particular, the scheme is very robust against the inefficiency of the 
detectors, which is clarified by a concise formula for the fidelity of the fixed point of 
the mixing map to the target singlet state. 

We here concentrated on a specific physical model with two qubits fixed along a 
ID channel, where ancilla qubits flow. The present analysis however provides a general 
guideline for turning a probabilistic convergence scheme into a scheme that works with 
probability arbitrarily close 1. In general, the former probabilistic scheme keeps the 
target state, as long as the measurement on an ancilla reports desired result. If the 
measurement outcome is not the desired one and if in such a case we are sure that the 
system is in an orthogonal state to the target, it is quite easy to design the feedback: 
we simply "shake" the system to provoke the transition from the orthogonal state to 
the target state. This inflow to the target state, in the absence of the outflow, ensures 
the convergence to the target state. Note also that the methods given in Ref. [TU] 
can be interpreted as feedback schemes for Ref. [2] . The methods developed for these 
simple systems pave the way for general strategies for driving systems to target states. 
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